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Abstract. We present experimental results showing the diffuse reflection of a 
Bose-Einstein condensate from a rough mirror, consisting of a dielectric substrate 
supporting a blue-detuned evanescent wave. The scattering is anisotropic, more 
pronounced in the direction of the surface propagation of the evanescent wave. 
These results agree very well with theoretical predictions. 
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1. Introduction 

The study of the interactions between ultra cold atoms and surfaces is of major interest 
in the context of Bose-Einstein condensation on microchips One motivation is to 

understand the limitations on integrated matter wave devices due to imperfect surface 
fabrication or finite temperature. For example, it has been shown that the quality 
of the wires used in microfabricated chips is directly linked to the fragmentation 
effects observed in Bose-Einstein condensates (BECs) trapped near a metallic wire 
. Moreover, the thermal fluctuations of the current in a metallic surface induce spin 
flip losses in an atomic cloud when the distance to the surface is smaller than 10 /im 
typically [1]. 

Dielectric surfaces and evanescent waves have also been explored for producing 
strong confinement. They have the advantage of a strong suppression of the spin 
flip loss mechanism compared to metallic structures jjj. With such a system, one can 
realize mirrors [HJ , diffraction gratings |Zj , 2D traps jHj or waveguides P; . Experiments 
involving ultra cold atoms from a BEC at the vicinity of a dielectric surface have 
recently made significant progress, leading for instance to the realization of a two 
dimensional BEC ^UJ) to the study of atom-surface reflection in the quantum regime 

and to sensitive measurements of adsorbate- induced surface polarization ^2] and 
of the Van der Waals/Casimir-Polder surface interaction 




Figure 1. Dielectric prism supporting the evanescent wave. The surface is coated 
by two layers of successively low and high refraction index to realize a wave guide 
and enhance the evanescent field. For each incident polarization, TE or TM, 
coupling is resonant for a given incident angle. The experiments were performed 
with TE polarization. One denotes x as the propagation axis of the evanescent 
wave along the surface, y as the other horizontal axis and z as the vertical one. 



In this paper, we present experimental results and a related theoretical analysis 
of Bose condensed Rubidium atoms interacting with the light field of an evanescent 
wave above a dielectric slab. The evanescent wave is detuned to the blue of an atomic 
transition line and provides a mirror for a BEC that is released from a trap and falls 
freely in the gravity field of the earth. After the bounce off the mirror, we observe 
a strong scattering of the atomic cloud (diffuse mirror reflection) that is due to the 
roughness of the slab surface where the evanescent wave is formed In our case, the 
phase front of the reflected matter waves is significantly distorted because the effective 
corrugation of the mirror is comparable to AdB/47rcos0 where AdB is the incident de 
Broglie wavelength and 9 the angle of incidence. This is similar to early experiments 
with evanescent waves |14| and with magnetic mirrors We mention that later 

experiments achived a significantly reduced diffuse reflection (Arnold et al and 
were even able to distinguish a specularly reflected matter wave (Savalli et al |17|). 
The key result of our experiment is that we can quantitivcly confirm the theoretical 
analysis developed by Hcnkcl et al 18 , combining independent measurements of the 
dielectric surface and the bouncing atoms. 

The paper starts with a presentation of the experiment and an analysis of the 
experimental results, following Ref.JIH|. We then outline an improved theoretical 
analysis based on Ref. JT8J and discuss the momentum distribution of the reflected 
atoms, in particular its diffuse spread and its isotropy. 



2. Setup 

The evanescent wave is produced by total internal reflection of a Gaussian laser beam 
at the surface of a dielectric prism. As shown in Fig^ the surface of the prism is 
coated by two dielectric layers, a TiC>2 layer on top of a Si02 spacer layer. This coating 
forms an optical waveguide that resonantly enhances the evanescent field above the top 
layer 120] ; we have designed this configuration for the study of two-dimensional atom 
traps [2] • The incident angle of the laser beam is fixed by the resonance condition for 
a waveguide mode; for the transverse electric (or s) polarization we use, the incident 
angle is di — 46.1° (at the TiO^/vacuum interface, index 7iTio 2 = 1-86). The resulting 
exponential decay length of the light field is = 93.8 nm, and I = Iq e~ 2KZ is the 
light intensity. 

The mirror light is produced by a laser diode of power 40 mW detuned 1.5 GHz 
above the atomic D2 line (A = 780 nm or 1/A = kj^jl-K = 12820 cm -1 ). The Gaussian 
beam is elliptical and produces on the surface a spot with l/\/e waist diameters of 
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220 film along x and 85 fim along y (see coordinate axes in FigQ. A measurement 
of the reflection threshold for the atom beam, taking into account the van der Waals 
attraction toward the surface (Landragin et al. in Ref.jHj), gives access to the light 
intensity at the surface in the spot center: Iq = 210 W/cm 2 . This value is lower than 
expected from the design of the dielectric coating; we attribute this to the losses due 
to the roughness of the deposited TiC>2 layer (see Figure |3] below and the discussion 
there) . 

3. Atom bounce 

3.1. Data 

The experiment proceeds as follows: approximately 10 s atoms are confined in the 
hyperfine ground state _F = 2,TO^ = 2ina Ioffe Pritchard (IP) type magnetic trap, 
3.6 mm above the evanescent mirror |21) . The magnetic trap is cigar shaped, x 
being its long axis. Oscillation frequencies are, respectively, lo x /2tt = 21 Hz and 
uj±/2it = 220 Hz in the radial directions (y and z). The atoms are evaporatively 
cooled to below the condensation threshold and about N = 3 x 10 5 atoms are released 
at t = by switching off the magnetic trapping fields. These atoms reach the mirror 
after free fall at t rc b = 27 ms and bounce on it with a velocity v t — 265mm/s (normal 
incidence 9 — 0, de Broglie wavelength AdB = lirh/mvi = 17.3 nm). Around the 
bouncing time i ra b, the mirror laser is switched on for At = 2.2 ms. Limiting this 
time window At prevents near-resonant photon scattering during free fall or after 
reflection. 

The atoms are detected by absorption imaging either before or after reflection. 
During free fall, the cloud expands along the radial directions because potential and 
interaction energy is released, but its width along x remains nearly constant. The 
analysis of pictures taken before reflection gives access to the following parameters: 
fraction of condensed atoms Nq/N = 0.4, kinetic temperature of thermal cloud 
T = 285 nK, initial Thomas-Fermi size along x of the condensed fraction R x = 90 fim 
and Thomas-Fermi velocity width along z: V± — 5.96 mm/s. The condensate velocity 
width along x is very small, thus non directly measurable. However, it can be inferred 
from the knowledge of V± and the oscillation frequencies in the magnetic trap, using 
the solution for an expanding BEC gU; we get = f ^-V± = 0.89 mm/s. The 
observation of the center of mass motion during free fall permits us to calibrate the 
pixel size knowing gravity's acceleration and to infer the initial position and velocity of 
the cloud. The magnetic field switching process communicates a small acceleration to 
the atoms along x, resulting in a horizontal velocity v x = —30.7 mm/s (see Figure|2J). 

After reflection, the absorption images change dramatically (figure^- The atoms 
occupy the surface of a scattering sphere, hence an elastic, but strongly diffuse 
scattering occurs. For t > t rc b, the cloud width along x increases from its initial 
value due to an additional velocity spread a Vlc . The velocity Gaussian radius at 1/ \fe. 
deduced from the pictures is 39.4 mm/s. Taking into account the initial velocity 
width before reflection, the spread due to diffuse reflection is a V;c — 39 mm/s, that 
is 6.6 ± 0.2 v TCC where u TO c = fi-fci/m = 5.89 mm/s is the recoil velocity for Rb. This 
corresponds to an angular (rms) spread A9 « 8.4°. 

The effect of diffuse reflection along y is more subtle to analyze, as this axis is 
aligned with the direction of observation. However, it is possible to extract information 
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Figure 2. Absorption imaging pictures of a bouncing BEC with 3 X 10 s atoms 
for different times of flight after reflection: 2 ms, 7 ms, 12 ms, 17 ms and 22 ms. 
The pictures, taken with a 200 /is long, resonant pulse, are merely superimposed. 
The wide black line in the bottom is due to the prism surface, slightly tilted from 
the imaging axis. Picture dimensions are 5.7 mm X 4.4 mm. 

about o~ Vy from the picture. If for instance the scattering were totally isotropic, with 
a v = o~ Vml the atomic cloud should extend asymmetrically towards — z at a given 
position x, as the projection of a spherical shell onto a plane extends towards the 
inner part of the circle (see figure 0J). If on the contrary the scattering would take 
place only along x, the cloud width along z at a given position x should be very small, 
with a symmetric shape. 

3.2. Simulation 

To get some insight into what happens along y, we performed a numerical simulation of 
the atomic reflection. The simulation calculates N — 3x 10 5 individual classical atomic 
trajectories. The initial positions and velocities are chosen to mimic the experimentally 
measured parameters: 40% of the atoms are "condensed" and are described by the 
initial 3D Thomas-Fermi velocity and position distribution. (We neglect the position 
spread along y and z because its contribution to the cloud size after a few ms of time 
of flight is very small.) The remaining 60% of the atoms are distributed according to 
gaussian profiles for velocity and position, with widths inferred from the knowledge of 
temperature and trap parameters. Position and velocity of the cloud centre are fixed 
to the experimental values as well. The mirror is modelled as an instantaneous diffuse 
reflector. This assumption is reasonable as the typical time spent in the evanescent 
wave is small, 1/nvi = 0.35 /is. After reflection, the atomic velocity is modified 
to describe both specular reflection (inversion of vertical velocity) and scattering. 
A random horizontal velocity is added to the reflected velocity with a gaussian 
distribution. We take a l/v^ radius o~ Vm = 39mm/s, as measured experimentally, 
and run simulations with varying a Vy . The z component of the velocity is adjusted 
in order to preserve kinetic energy (the scattering process is elastic, total energy is 
conserved). The simulation also takes into account spontaneous emission. For our 
parameters, the atom spontaneously emits on average 0.13 photons per bounce |23|. 
We randomly draw the number of photons from a Poisson distribution and add a 
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Figure 3. Simulation of a bouncing BEC with 3 X 10 5 atoms for the same times 
of flight than the experimental ones, figure [5] For this series, the velocity spread 
was chosen to be a Vy = (t Vx /2 = 19.5 mm/s. The position of the mirror surface 
is marked by a grey line. 



recoil of 1 v rcc in a random direction in velocity space for each emission event. After 
calculation of all atomic trajectories, the atomic density profile is integrated along y as 
in the experimental pictures. We finally apply a Gaussian blur filter (width <7 ros = 9 /im 
along x and 20 fiiii along z) to mimic the finite resolution of the experimental imaging 
system that we calibrated independently. 

3.3. Anisotropic scattering 

The qualitative agreement between the experimental and simulated pictures is very 
good as can be seen on figure [3 To be more quantitative for the possible values of 
the velocity spread a v , we analyze the central part of the cloud. For each time of 
flight, a region of size 0.8 mm x 1.5 mm along x and z respectively, centred on the 
maximum density of the cloud and identical for experimental and simulated pictures, 
is isolated and an integration of the signal is performed along x. We are left with a 
cut of the cloud along z, averaged over 0.8 mm along x. The experimental profile is 
compared to the simulated one, for different choices of o~ Vy after the bounce. Results 
are shown on figure 0] for a time of flight 59 ms. 

The experimental data clearly exclude an isotropic diffuse reflection (figure 0J 
bold line) . They also are different from the pure one dimensional scattering case (thin 
line): what fits best of all is a model intermediate between these two extremes, i.e. the 
scattering is only half as strong along y compared to x. The atom mirror thus has an 
angular reflection characteristic that is elongated in the direction parallel to the (real 
part of the) wave vector of the evanescent wave. Spontaneous emission plays only a 
minor role for our parameters, but we found that the agreement with the experimental 
density profiles is improved by taking it into account, in particular on the lower left 
wing of the peak. 
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Figure 4. Atomic density profiles integrated along y and averaged along x, after 
59 ms total time of flight, i.e. 32 ms after reflection. Closed circles: normalized 
experimental data. Lines: result of numerical calculation with a starting height 
3.59 mm above the mirror, N /N = 0.4, V v = V z = V± = 5.96 mm/s, 
V x = 0.89 mm/s, R x = 90 /im, T = 285 nK, v x = -30.67 mm/s, u z = 0.3 mm/s 
and a Vn . = 39 mm/s, values deduced from the experimental pictures. Thin line: 
totally anisotropic scattering (& Vy = 0); dashed line: anisotropic scattering 
with a Vy = <r„ x /2; bold line: isotropic scattering (cr Vy = cr V:c ). All curves are 
normalized to a maximum value of unity. 



3.4- Mirror corrugation 

For a theoretical prediction of the anisotropic mirror reflection, we use the theory 
of Ref.^21 where the diffuse scattering is attributed to the interference between the 
evanescent wave and light diffusely scattered from the rough glass surface. Within 
this theory, one can compute the width of the momentum distribution of the reflected 
atoms provided the power spectrum of the surface roughness is known. This power 
spectrum is a quantitative measure of the surface quality and has been measured with 
an atomic force microscope (AFM). A typical 4.5 x 4.5^m 2 portion of the surface of 
the coated prism is shown in figure One sees the top face of pillar-like structures 
which are typical for epitaxially grown TiC>2 on a substrate. The AFM data yield a 
surface roughness a — 3.34 nm (the rms spread of the measured surface profile). A 
Fourier transform of the AFM image gives access to the power spectrum Pg(Q). (We 
use capitalized boldface letters for two-dimensional vectors in the mirror plane.) ft is 
found to be isotropic (a function of Q only) and well fitted in the wave vector range 
1 . . . 13 Ul by a power law with a low-frequency cut-off (see figure 

The fit gives access to the parameters a — 4.8, Pq = 5.3 x 10~ 4 k L A and Qo = 4.94 Ul- 
In terms of this power spectrum, the rms surface roughness a is given by 
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Figure 5. (left) Typical AFM picture of the prism surface. The dimensions are 
4.5 fimx4.5 fim. The grains are the top facets of pillar-like structures characteristic 
of an epitaxed Ti02 surface. 



Figure 6. (right) Open circles: Power spectrum Ps(Q) deduced from the AFM 
measurement. Solid line: fit of the data by the function given at Eq.ITt, 



and the fitted parameters yield a = 3.36 nm, in excellent agreement with the value 
directly deduced from the rms spread of the AFM data. 

3.5. Diffuse reflection theory and comparison to the data 

We now show that the diffuse reflection we observe can be understood within the 
theory of Ref. ^H] . We note first that the cloud is so dilute at the bounce that a single- 
atom picture is sufficient to capture the physics |24| . For a fixed incident momentum 
Pi nc = P — e z mvi near normal incidence, the reflected wave function can be written 
in the form of a plane wave with a randomly modulated phase front: 



where N is a normalization factor and p spe c = P + e z mvi. (The dependence on the 
angle of incidence is actually negligible for our parameters US]-) The phase 5c/)(R) 
depends on the 'impact position' R on the mirror, i.e., the projection of r onto the 
mirror plane. We perform an ensemble average over the realizations of the rough 
surface and compute the atomic momentum distribution Pa(P +hQ) from the (spatial) 
Fourier transform of the 'atomic coherence function' (Sec. 6 of Ref.|18j1 



(We take (<5</>(R)) = 0, assuming the roughness to be statistically homogeneous.) The 
variance of the phase shift can be found from the following formula (Eqs.(6.15) and 



^refl(r) = iVexpi [p spoc • r + 8<f)(R)} 



(3) 



(V* efl (r)Vrefl(r')) = ^ 2 expi [p spcc • (r' - r)] 

xexp -i((ty(R)-ty(R')) 2 ) 



(4) 
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(5.16) of Ref.dSj) 




(5) 



where -B a t(Q) is the "atomic response function" given in Eq.(5.15) of Ref.|18|. For 
the parameters of our experiment, we find that the phase shift has a variance 
((50 2 (R)) = 16.5 large compared to unity. In this regime, Ref.^2] nas shown that 
the reflected atomic velocity distribution approaches a Gaussian shape whose width 
along the cc-direction, for example, is given by 



This expression gives the additional broadening of the incident velocity distribution 
due to the diffuse mirror reflection. We perform the integration of Eq.@ numerically, 
with the roughness power spectrum determined previously from the AFM images 
(Eq.J3>). For simplicity, we calculate the response function -B a t(Q) using scalar light 
scattering from the topmost interface only, ignoring the actual layered structure. We 
believe that this approximation is sufficient, at least for describing the scattering in 
the x-direction: as shown in Ref |18| . the atom does not change its magnetic sublevel 
if it scatters in this direction and if the evanescent wave is linearly polarized. These 
conditions are met here so that both atom and light can be described by scalar wave 



Within the theoretical model outlined above, the velocity spread along the 
propagation direction of the evanescent wave is found to be <t„ x = 6.76 v rcc . This 
value is in very good agreement with the experimental value 6.6 ± 0.2 v rcc . This is 
a very satisfying result because the theory only contains, within the approximations 
we made, parameters that are based on independent measurements. We believe that 
this is the first quantitative demonstration of evanescent wave scattering in the diffuse 
regime. 

3.6. Discussion of the anisotropy 

We also compute the anisotropy of the reflected atoms and find a ratio cr Va ./a v = 2.6, 
in good agreement with the value (2 ± 0.5) extracted from the experimental data. As 
discussed in Ref. ^Hl i this anisotropy arises from the fact that diffuse reflection occurs 
predominantly by Bragg transitions where a photon is absorbed from the evanescent 
wave (with wave vector k x = fcL"-Ti0 2 sin ) and another photon is emitted into a 
diffusely scattered mode that emerges at grazing incidence into the vacuum half- 
space (or the inverse process). If these scattered modes are distributed isotropically 
in the mirror plane on a circle of radius r sc fci, the ratio of the rms spreads would 
be a v ^ja Vy = (2(nxi0 2 sin(#j)/r sc ) 2 + l) 1 / 2 . Taking r sc = 1, which corresponds to 
scattered modes emerging at grazing incidence, we again find an anisotropy ratio 
of w 2.5. This agreement is not very surprising since the rough surface has a 
power spectrum much broader than the photon wavenumber (FiglHJl. Within this 
simple calculation, however, we can also get a quick estimate of the impact of the 
dielectric coating. The choice r sc = riTi0 2 sin(0j) corresponds to resonant scattering 
into waveguide modes in the TiC>2 layer and leads to a ratio o~ V:c /a Vy ~ a/3 which 
cannot be excluded experimentally. 




(6) 



fields. 
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In conclusion, we have observed the diffuse reflection of an ultracold atomic beam 
from an evanescent wave. The wave propagates on the rough surface of a dielectric 
prism, and light scattering leads to an atom mirror showing a significantly nonspecular 
reflection. The angular broadening of the reflected atoms, as well as their anisotropic 
angular distribution in the mirror plane, are in good agreement with a theory 
developed by Henkel et al. It is remarkable that this agreement does not 

imply any free parameters since we independently measured the spectrum of the 
surface roughness with an AFM. In our experiment, using a BEC has mainly practical 
advantages. Indeed, as we mentioned above, everything can be understood within a 
single-atom picture, and after diffuse scattering, spatial coherence is seriously reduced, 
as is discussed in Ref.^B] and investigated in Ref.jjS]. Nevertheless, the BEC provides 
crucial advantages because we achieve a very clean situation. Apart from a very low 
velocity spread V x <§C a Vw , a BEC has a negligible size when impacting the evanescent 
wave surface. This removes the need to take into account the mirror curvature due to 
the gaussian spot profile; the contribution of the initial size to the cloud width after 
reflection is negligible; and the losses given the finite size of the mirror (the waist of 
the reflected laser beam) are minimal. In fact, with a freely falling, ultracold, but 
thermal gas, the finite mirror size would lead to strongly reduced signal. 
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